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ABSTRACT 



Fibromyalgia syndrome is a chronic pain disorder characterised by widespread pain and tenderness in muscles 
and deep tissues. Current theories regarding the pathophysiological origins of fibromyalgia syndrome point 
towards central sensitisation and a decreased capacity of descending nociceptive controls. Morphological 
alterations to subcortical brain regions may contribute to such pathophysiological mechanisms, and to pain 
and other symptoms seen in fibromyalgia. Therefore, we evaluated geometric differences in subcortical 
structures in fibromyalgia patients relative to healthy people using a novel method of shape analysis. Sixteen 
female fibromyalgia patients and 15 age and sex matched, healthy control subjects underwent high- 
resolution Tl -weighted magnetic resonance image scanning. Data was analysed using shape analysis of 
15 subcortical regions and standard voxel-based morphometry analysis. 

Fibromyalgia syndrome patients, relative to healthy control participants, exhibited alterations to the shape of 
the left lateral aspect of the lower brainstem (medulla). The mean total volume of the brainstem was also 
found to be significantly reduced in the patient group compared to healthy control subjects, and this brainstem 
volume reduction in patient group significantly correlated with clinical manual tender point scale scores. Voxel- 
based morphometry analysis revealed that patients also demonstrated decreased local grey matter volumes in 
the brainstem (pons) and left precuneus, and increased grey matter volumes in bilateral primary somatosensory 
cortices. 

Results suggest that the volume reduction and associated geometric shape alterations seen in the brainstem of 
the patient group may contribute to sensitivity to pressure pain in fibromyalgia syndrome. This finding may be 
due to structure-related deficiencies in regions subserving descending nociceptive control. 

© 2013 The Authors. Published by Elsevier Inc. All rights reserved. 



1. Introduction 

Fibromyalgia syndrome (FMS) is a chronic pain disorder characterised 
by widespread pain and tenderness, morning stiffness, sleep distur- 
bance and fatigue (Bennett et al., 2007; Wolfe et al, 1 990). FMS patients 
also demonstrate hyperalgesia and allodynia and cognitive deficits 
(Clauw, 2009) and dysfunctional postural control (Jones et al, 2011). 
The disorder is often co-morbid with other diagnoses such as chronic 
fatigue syndrome, temporomandibular joint disorder and irritable 
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bowel syndrome (Clauw, 2009). In 1990, the American College of 
Rheumatology (ACR) empirically defined FMS as pain during manual 
palpation of at least 1 1 out of 18 designated tender points and affecting 
all 4 quadrants of the body for a period of at least three months (Wolfe 
et al., 1990). A more recent ACR criterion, without a need for a physical 
examination, classifies FMS patients according to widespread pain 
index and symptom severity scores (Wolfe et al, 2010). 

As alterations to central processing of somatosensory stimuli may 
influence symptoms of FMS (Clauw et al, 201 1 ), studies have previously 
utilised voxel based morphometry (VBM) to investigate structural grey 
matter differences in FMS patients relative to healthy control subjects 
(Burgmer et al, 2009; Kuchinad et al, 2007; Luerding et al, 2008; 
Lutz et al, 2008; Robinson et al, 2011; Schmidt-Wilcke et al, 2007; 
Valet et al, 2009). FMS patients have demonstrated local grey matter 
density reductions in cingulate, insular cortices, medial frontal cortices, 
parahippocampal gyri (Kuchinad et al, 2007), superior temporal gyrus 
and left posterior thalamus (Schmidt-Wilcke et al, 2007). Grey matter 
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density increases were shown in the left orbitofrontal cortex, cerebel- 
lum and bilateral striatum ( Schmidt- Wilcke et al, 2007). Analyses 
excluding age and depression identified reduced grey matter volumes 
in prefrontal, cingulate and insular cortices of FMS patients (Valet 
et al, 2009), although an alternative study suggested that alterations 
were negated when implementing depression as a covariant (Hsu 
et al, 2009). VBM analyses of pre-defined regions of interest of FMS 
patients have shown reduced grey matter volumes in the cingulate, pre- 
frontal and insular cortices, hippocampi and amygdala (Burgmer et al., 
2009; Lutz et al., 2008; Robinson et al., 2011). 

Unfortunately, there is little consistency in the VBM data from FMS 
studies. This may result from methodological issues pertaining to tissue 
classification and arbitrary smoothing (Jones et al., 2005; Patenaude 
et al., 2011; Smith et al., 2006). Factors such as age, sample size and 
relative proportions of symptoms contributing to the complex symptom 
profile in the FMS population can also affect findings (May, 201 1 ). VBM 
has also been used to identify structural differences in patients with 
irritable bowel syndrome (Seminowicz et al., 2010), tension type 
headache (Schmidt- Wilcke et al, 2005), chronic fatigue syndrome 
(de Lange et al, 2005) and post-traumatic stress disorder (Chen 
et al., 2006; Villarreal et al., 2002). FMS is often comorbid with 
these disorders, and the overlap of symptoms has led some researchers 
to propose that similar neurophysiological mechanisms may be 
involved (Ablin and Clauw, 2009; Phillips and Clauw, 2011; Williams 
and Clauw, 2009). It is possible that a particularly high prevalence of a 
specific comorbidity in a heterogeneous FMS population could 
drive regional grey matter changes identified using VBM and explain 
the wide variety in the data. 

In the present study, shape analysis of fifteen subcortical regions was 
carried out in a homogenous sample of female FMS patients and 
healthy, age-matched control subjects. This method can be used to com- 
paratively investigate shape and volume of subcortical structures, and 
has recently been used to investigate abnormalities in long term absti- 
nent alcoholics (Sameti et al., 2011), patients with Alzheimer's disease 
(Zarei et al, 2010) and age-related changes in healthy populations 
(Goodro et al., 2012). This novel technique has been recently validated 
and shown to be consistent over a variety of magnetic field strengths 
and acquisition systems (Goodro et al., 2012). Shape analysis can locate 
subcortical structural boundaries without the need for smoothing. 
Therefore, it can accurately identify both location and direction of com- 
plex morphological alterations through direct measurement of geomet- 
ric shape more precisely than VBM (Patenaude et al., 2011). In order 
to compare results with shape analyses concurrent global and local 
grey matter alterations throughout the whole brain of FMS patients 
were also investigated using VBM. We hypothesised that FMS patients 
would show subcortical abnormalities in shape and volume and that 
the degree of such morphological changes would correlate with clinical 
measures of FMS severity such as manual tender point scale (MTPS) 
scores (Wolfe et al., 1990), Beck Depression Inventory (BDI) scores 
(Beck et al., 1961) and Fibromyalgia Impact Questionnaire (FIQ) scores 
(Burckhardt et al., 1991 ), which included self-report measures of dura- 
tion of symptoms. 

2. Methods 

2.1. Patients 

Sixteen female patients (age 38.5 ± 8.45 years, mean ± SD) took 
part in the study. Mean duration of symptoms was 9.1 ± 6.8 years, 
and mean time since diagnosis was 2.8 ±1.3 years (mean ± SD). All 
patients fulfilled ACR criteria for diagnosis with fibromyalgia on the 
day of scanning (Wolfe et al., 1990). Patients were recruited from 
two outpatient fibromyalgia clinics at regional NHS Foundation 
Trust hospitals; the Walton Centre, Liverpool, United Kingdom, and 
Wirral University Teaching Hospital, Wirral, United Kingdom. Informed 
consent was obtained from all participants in accordance with the 



Declaration of Helsinki and the study was approved by the National 
Research Ethics Committee of the United Kingdom, and the Research 
Governance Committees of both hospitals. Patients aged between 22 
and 52 years were considered for participation. Those with additional 
disease or disorders (not commonly comorbid with FMS such as diabe- 
tes or hypertension) were excluded, as were patients with a past history 
of major disease, alcohol/drug abuse or serious head or brain injury. 
Patients using analgesics (such as co-codamol) withdrew for 3 
to 5 days prior to scanning. For example, a patient taking 6 to 8 tablets 
of mild co-codamol (8/500 mg) a day was asked to discontinue 
use for 3 days, a patient using up to 8 tablets of strong co-codamol 
(30/500 mg), or high doses of dihydrocodeine, was asked to taper the 
dose over a period of 2 days before discontinuing use for 3 days prior 
to testing; withdrawal was managed by the clinical team during consul- 
tation. Analgesic medications with minimal central nervous system 
effects, such as paracetamol, were permitted. Patients using medica- 
tions with central nervous system effects (such as trycyclic antidepres- 
sants), which were not deemed suitable for withdrawal by the clinical 
team, were excluded. Common FMS medications at stable low doses 
such as pregabalin (up to 75 mg twice a day), gabapentin (up to 
300 mg twice a day) and amitriptyline (up to 10 mg at night), were 
considered acceptable for the criteria of minimal central nervous effects. 

2.2. Controls 

Fifteen age-matched female controls (age 39.4 ± 8.7 years, 
mean ± SD) were recruited through campus advertisement. Volun- 
teers taking regular medication, currently diagnosed with any disease 
or disorder, or demonstrating a history of major disease, alcohol/drug 
abuse or serious head or brain injury were excluded. All patients and 
volunteers were compensated for time and travel expenses. 

2.3. MRI data acquisition 

High-resolution 3-dimensional Tl -weighted images of the brain 
were acquired using a whole-body Trio 3 T MRI imaging system 
(Siemens, Erlangen Germany) and an 8-channel head coil. A Tl 
weighted modified driven equilibrium Fourier transform (MDEFT) 
sequence (TR = 7.92 ms, TE = 2.48 ms, flip angle = 16°, 176 sagittal 
slices, slice thickness 1 mm, matrix 256 x 256, in-plane voxel size 
1 mm x 1 mm, total acquisition time 12:51 min) was utilised. Fol- 
lowing the MR scan all participants underwent a clinical MTPS exam- 
ination (Wolfe et al., 1990), eighteen anatomically standard FMS tender 
points were palpated for 4 s using the thumb pad of the dominant hand. 
Pressure began at 1 kgf and was incremented by 1 kg per second until a 
maximum pressure of 4 kg is achieved. Following examination of each 
point patients reported whether they felt any pain during palpation 
and rated the pain verbally on a scale of 0 (no pain) to 10 (worst pain 
ever experienced). All participants also completed the BDI and FIQto 
evaluate affective disturbance and the impact of symptoms on quality 
of life in the week preceding the scan. Questionnaires and MTPS exam- 
inations were administered immediately following the MR scanning 
session (on the same day), so as to reflect medication levels and symp- 
toms appropriate to the time of scanning. 

2.4. Subcortical segmentation, vertex and volumetric analysis 

Fifteen subcortical structures (brainstem, bilateral thalami, hippo- 
campi, amygdalae, putamen, caudate nucleus, accumbens nucleus 
and pallidum) were segmented from each subject's high-resolution 
Tl structural scan using the Oxford Centre for Functional MRI of the 
Brain's (FMRIB) integrated registration and segmentation tool (FIRST) 
toolbox (Patenaude et al, 2011) in the FMRIB software library [FSL, 
http://www.fmrib.ox.ac.uk/fsl, (Smith et al., 2004)]. This method imple- 
ments a probabilistic adaptation of the active appearance model using 
a large, manually labelled data set as a training template wherein the 
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subcortical structures are parameterised as surface meshes with 
established vertices. In an automated process, the images are linearly 
registered (12 degrees of freedom) to the training template and subse- 
quently to a standard Montreal Neurological Institute (MNI) space sub- 
cortical mask template in FSL Finally, shape and intensity variations in 
the images are used to automatically segment structures in native 
space using the inverse transformation. 

To evaluate group differences in the shape of subcortical structures 
vertex analysis was performed. A multivariate Gaussian model of vertex 
location and intensity variation was used to generate a surface mesh for 
each structure in each subject. The number and correspondence of mesh 
vertices remain constant across subjects to allow for point-to-point 
group comparisons of mean vertex positions. Corresponding vertices 
of surface meshes for each individual subject were compared in stan- 
dard MNI space using the Multivariate General Linear Model with 
Pillai's trace as the test statistic. This results in a multivariate F-statistic 
for each vertex which is sensitive to between group differences in vertex 
coordinates in any direction (Patenaude et al, 2011). Geometric shape 
analysis of subcortical structures may be preferable to VBM for identifi- 
cation of subtle morphological alterations as this method identifies the 
location and direction of complex morphological alterations through 
direct measurement of shape without the need for smoothing 
(Patenaude et al., 2011). The accuracy of this method has been shown 
to be superior to previous automated segmentation tools (Patenaude 
et al, 201 1 ), and FIRST analysis is consistent over a variety of magnetic 
field strengths and acquisition systems (Goodro et al., 2012). This method 
has also been compared to the 'gold standard' of segmentation in struc- 
tural imaging, manual segmentation of images by an expert, and the find- 
ings indicate a high degree of accuracy (Babalola et al., 2009). Individual 
vertex results were corrected for multiple comparisons (across vertices) 
using false discovery rate (FDR) method (Benjamini and Hochberg, 
1995) and P < 0.05 was considered significant. Bonferroni-Sidak correc- 
tion was also employed to account for the multiple comparisons required 
across 15 structures. 

As recent studies pointed to an association between Chiari I malfor- 
mation and FMS (Heffez, 2011; Thimineur et al, 2002; Watson et al., 
2011), each anatomical image was assessed for Chiari I malformation 
by measuring the distance from the inferior tip of the cerebellar tonsils 
to the line connecting basion and opisthion (Watson et al., 2011). This 
measure was of importance for evaluation of potential shape or volume 
differences because increased pressure of cerebrovascular fluid below 
foramen magnum in Chiari I malformation may affect the shape of 
subcortical structures such as the brainstem. 

For structures demonstrating a significant shape alteration in FMS pa- 
tients relative to healthy participants, multivariate analysis of covariance 
(MANCOVA) was performed using clinical and psychological measures 
as covariates to investigate whether clinical or psychological variables 
would explain group shape differences. However, if the measures 
violated the assumption of homogeneity of regression slopes (i.e. dem- 
onstrated a bimodal distribution which was confirmed by a significant 
difference in group means using an independent samples t-test), then 
they were deemed unsuitable for use in MANCOVA analysis. Instead, 
these measures were implemented in a correlation analysis using FIRST 
in the FMS patient group only. These variables were independently 
entered as correlates in the design and vertex analysis was performed 
to evaluate whether they correlated with vertex positions in structures 
showing significant group differences in shape. FDR correction was 
applied to correct for multiple tests across vertices. Chiari I parameters, 
MTPS scores, FIQ scores, BDI scores and individual brain structure vol- 
umes were the measures utilised as potential confounds in MANCOVA 
or correlation analyses to investigate significant group differences. 

The surface meshes generated for each subject were then 
transformed back into native space and boundary corrected before 
the volume (mm 3 ) for each structure was calculated. For each subcorti- 
cal structure demonstrating a significant difference in shape analysis, 
volumetric data were exported for each participant and mean values 



for FMS patient group and healthy control group were compared using 
Student's independent sample f-tests. In the event of significant group 
differences in the volume of specific structures, univariate ANCOVA was 
performed implementing total intracranial volumes and total grey mat- 
ter volumes from VBM analyses as covariates to evaluate whether such 
factors would account for the findings. Individual volumes in the FMS 
patient group were also utilised for Pearson's correlation analyses with 
clinical measures (MTPS scores) in SPSS v.19 (SPSS Inc., Chicago, USA). 

2.5. Voxel based morphometry analysis 

VBM pre-processing and analysis were performed using VBM8 tool- 
box (http://dbm.neuro.uni-jena.de/vbm) in SPM8 (Welcome Trust Cen- 
tre for Neuroimaging, University College London, United Kingdom) 
in Matlab v.7.8 (The MathWorks Inc., USA). Voxel based morphometry 
consists of several steps including spatial normalisation of images, seg- 
mentation of grey matter and other tissues, and analysis of global and 
local differences in voxel intensities. Images were registered to MNI 
space and an optimized method of VBM using the VBM8 Toolbox was 
implemented for segmentation of grey and white matter, default set- 
tings were used unless otherwise indicated. A hidden Markov random 
field model is applied as part of this procedure to reduce noise effects 
and, as spatial normalisation causes expansion or contraction of some 
brain regions, segmented images were scaled by the degree of contrac- 
tion (modulated), so that the total volume of grey or white matter in the 
images remains the same as in the original images. Images were 
smoothed using an isotropic Gaussian kernel of 10 mm full width at 
half maximum. A data quality check based on inhomogeneity measures 
of the sample, as implemented in VBM8 toolbox, was used to check for 
anomalous data (outside 2 standard deviations) which was visually 
inspected and subsequently excluded if necessary. Using these criteria 
no data was excluded from the sample. 

A voxel-wise comparison was performed between the FMS patient 
and healthy control group using the general linear model implemented 
in VBM8. As in previous clinical studies (Buckalew et al, 2008; Schmidt- 
Wilcke et al., 2006) where pain-related grey matter alterations may be 
small, whole brain differences in voxel intensities were first evaluated 
using an initial threshold of P < 0.001 (uncorrected) to identify clusters 
of contiguous voxels (thresholded to encompass at least 30 voxels) 
showing a difference. Subsequently, the clusters of voxels demon- 
strating grey matter volume differences at the uncorrected level were 
compared using a statistical threshold corrected for family-wise error 
(FWE) at cluster level (P < 0.05 corrected). Finally, for any regions 
demonstrating significant group differences, individual regional grey 
matter volumes in the FMS patient group were exported and utilised 
for Pearson's correlation analysis with MTPS, FIQ and BDI scores in 
SPSS v.19. This analysis was performed to evaluate whether clinical or 
psychological factors would account for alterations to regional grey 
matter volumes in the FMS group. 

2.6. Total intracranial tissue volumes and correlations with clinical measures 

Total volumes of grey matter, white matter and total intracranial 
volumes for each participant were exported and mean values for 
FMS patient group and healthy control groups were compared using 
Student's independent sample f-tests. Pearson's correlation analysis 
was performed to evaluate possible correlations between total grey 
matter volumes in the patient group and clinical measures of symptom 
severity (MTPS, BDI and FIQ scores). 

3. Results 

3.1. Clinical and psychological analysis 

A Student's independent t-test indicated a significant difference be- 
tween the mean MTPS scores for patients (4.98 ± 1.89) and controls 
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(0.26 ± 0.34); t (29) = 9.48, P < 0.001. Mean FIQscores in the FMS pa- 
tient group (63.37 ± 15.83) also differed significantly from those in the 
healthy control participants (6.26 ± 6.51); t (29) = 12.74, P < 0.001. A 
significant difference was also seen between mean BDI scores for FMS 
patients (19.50 ± 11.19) and controls (3.87 ± 4.42); t (29) = 5.05, 
P< 0.001. 

32. Shape analysis 

Vertex analysis was performed to evaluate the shape and volume of 
1 5 subcortical structures in FMS and healthy control groups. However, 
only the brainstem showed a localised shape difference, which is signif- 
icant following FDR correction for the number of vertices tested and 
Bonferroni-Sidak correction for multiple comparisons across structures 
(P = 0.01 corr)- In the FMS patient group the brainstem demonstrated an 
inward movement of vertices compared to healthy control participants 
on the left lateral medullary funiculus extending from the inferior tip of 
the brainstem to the level of the inferior olivary body (Naidich and 
Duvernoy, 2009) suggesting a volume reduction in this region. A signif- 
icant outward movement of fewer vertices located above this area in the 
left medulla at the level of the inferior olivary body was also apparent 
which further suggests an alteration to geometric shape in the region. 
Taken together these results suggest an inward shape change to the 
left lateral aspect of the brainstem in FMS patients at the level of the me- 
dulla accompanied by a smaller (in size) outward shape change at the 
level of the inferior olivary body. Fig. 1 shows the locations, directions 




Fig. 1. Vertex analysis of shape alterations in the brainstem of the FMS patient group in 
comparison to healthy control participants. Upper panel shows the anatomical location 
of the brainstem and the local area exhibiting shape change in FMS patient group. Lower 
left panel indicates shape change in FMS patient group compared to healthy control 
group at an uncorrected level, this transparent image shows direction of vectors, inward 
direction represents relative inward positions of vertices (in FMS subjects compared 
with healthy control subjects) indicative of volume reduction. Outward direction of 
vectors indicates relative shape increases in FMS patient group. Arrow colour and surface 
colour indicate the F-statistic of the change in the specific vertices (see colour bar). Lower 
centre panel shows the location of the difference in patient group following FDR correction, 
red colour indicates areas which did not differ significantly from healthy control group. 
Lower right panel shows a semi-transparent image following FDR correction, indicating 
the direction of significant vectors showing alterations in the FMS group in comparison 
to healthy control participants. (For interpretation of the references to colour in this figure 
legend, the reader is referred to the web version of this article.) 



and F-statistics of shape change of brainstem vertices in the FMS 
group relative to healthy control group. 

A Student's independent t-test was performed to compare mean 
brainstem volumes, as indicated by FIRST analysis, in FMS patient and 
healthy control groups. The mean total volume of the brainstem in 
FMS patient group was shown to be significantly reduced in comparison 
to healthy control participants (t (29) = 2.56, P = 0.016). Fig. 2A 
shows the bar charts and error bars for both group's mean brainstem 
volumes. Univariate ANCOVA analysis revealed that the group differ- 
ence in volume of the brainstem was similarly significant following 
correction for total intracranial volume and total grey matter volume 
covariates (F(l,27) = 7.30, P = 0.012). Pearson's correlation analysis 
was performed to evaluate the relationship between the reduction of 
brainstem volume and MTPS scores, BDI scores and FIQ scores in the 
FMS patient group. A significant one- tailed correlation (r = —0.45, 
N = 16, P = 0.039) was evident, indicating that patients exhibiting 
greater reductions in volume of the brainstem reported higher scores 
for the MTPS evaluation. Fig. 2B shows the scatter plot of individual 
brainstem volume and MTPS score data for the FMS patient group dis- 
tributed along the regression line. 

As Chiari I malformation might contribute to the shape and volume 
changes in the brainstem, we analysed the distance between the 
basion-opisthion line and lower tip of cerebellar tonsils (Watson 
et al., 201 1 ). Two FMS patients and one control subject showed a tonsil- 
lar position >5 mm below the basion-opisthion line which meets the 



A 

Brainstem 
[mm3] 

24000 




20000 
19000 
18000 



B 

Brainstem 
[mm3] 

24000 



22000 



20000- 



18000 



16000 



Controls 



Patients 




MTPS 



Fig. 2. A: Bar chart illustrating the mean total volume of the brainstem in FMS patient 
group and healthy control participants and error bars. * indicates a significant difference 
in an independent Student's t-test comparison of group mean values. B: Scatter plot 
showing manual tender point scale (MTPS) scores and total volume of the brainstem 
(mm 3 ) as calculated during FIRST analyses. The linear regression line is also shown. 
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diagnostic criterion of Chiari I malformation. The mean distance be- 
tween the basion-opisthion line to the tip of cerebellar tonsils was 
0.58 ± 3.33 mm (mean ± SD) in FMS group and -0.12 ± 2.23 mm 
(mean ± SD) in healthy control subjects (t (29) = -0.68, P > 0.05), 
positive values denoting measurements inferior to the basion-opisthion 
line. This distance measurement was entered as a covariate in a 
MANCOVA design vertex analysis in FIRST. The vertices differentiating 
FMS patients and controls in the lower brainstem (Fig. 1A) remained 
unchanged after entering the distance between basion-opisthion line 
and the lower tip of cerebellar tonsils as a covariate. Thus, it is unlikely 
that Chiari I malformation would account for the reduction of lower 
brainstem volume seen in FMS patients. The additional clinical and psy- 
chological data collected during the study demonstrated a significant 
bimodal distribution and was therefore not appropriate for MANCOVA 
analysis. These variables (MTPS scores, BDI scores, FIQ scores and indi- 
vidual brainstem volumes) were independently employed in correla- 
tion analysis using FIRST in the FMS patient group. No significant 
correlations were identified between vertex positions in the brainstem 
and any of the clinical or psychological measures in the FMS patient 
group. 

3.3. Regional grey matter volumes 

Voxel-wise comparison of local grey matter volumes across the 
whole brain revealed that FMS patients, relative to healthy control par- 
ticipants, exhibited two significant clusters of grey matter volume 
reduction. One was located in the brainstem in the left ventral aspect 
of the basilar pons, in the pontine nuclei, and the other in the left 
precuneus (cluster extent P < 0.05, FWE-corrected). FMS patients also 
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Fig. 3. A: Local grey matter volume decreases as indicated by clusters representing spatial- 
ly extended groups of voxels which differed significantly in grey matter volume at the 
uncorrected level (P < 0.001) in whole brain analysis and corrected level (P < 0.05) in 
cluster level analysis. Clusters are displayed in glass brains {upper panel) and MNI 
standardised anatomical brains (lower panel) x, y and z co-ordinates indicate slice dimen- 
sions in MNI space. B: Local grey matter volume increases at the same significance level. 
Clusters are displayed in glass brains (upper panel) and MNI standardised anatomical 
brains (lower panel). 



Table 1 

Local grey matter volume decreases in FMS patient group relative to healthy control group. 



Structure 




MNI [mm] 




k 


Z 


T 


Sig. 




x 


y 


z 










Brainstem 


-14 


-21 - 


-38 


35 


3.41 


3.82 


0.008 


Precuneus 


-23 


-51 


12 


40 


3.37 


3.77 


0.008 



exhibited two clusters of grey matter volume increases located in 
bilateral primary somatosensory cortices. Fig. 3 shows the locations of 
clusters of voxels demonstrating grey matter alterations in FMS patient 
group relative to healthy control subjects displayed in glass brains 
and MNI standardised anatomical brains. Table 1 shows the MNI co- 
ordinates of clusters demonstrating grey matter volume decreases. 
Table 2 shows increases in FMS patients relative to healthy control sub- 
jects. The anatomical locations of t-maxima as defined by the Harvard- 
Oxford atlas in FSL, peak T and Z values, number of voxels (k) and 
the cluster-level FWE corrected significance values are also shown. 
Pearson's correlation analysis was used to investigate potential linear 
relationships between regional grey matter volume changes and 
MTPS, BDI and FIQ scores within the FMS patient group. Grey matter 
volume data was extracted for all clusters demonstrating increases or 
decreases for each patient in the FMS group and Pearson's correlations 
with MTPS scores, BDI scores and FIQ scores were calculated. No signif- 
icant correlations were found. 

3.4. Total intracranial volume statistics and correlations with clinical 
measures 

A Student's independent samples t-test indicated no significant dif- 
ferences in mean total grey matter volume for FMS patients compared 
to healthy control subjects (P > 0.05), although total mean grey matter 
volume in FMS patient group (626.30 cm 3 ) was lower than in the healthy 
control group (644.19 cm 3 ). Pearson's correlation analysis showed a 
significant correlation between total grey matter volume and MTPS 
scores in the FMS patient group (r = -0.63, N = 16, P = 0.009). 
Patients with lower total grey matter volumes scored higher on the 
MTPS examination. Fig. 4 shows the data distributed along the regression 
line in a linear relationship. No correlations were seen between total 
intracranial tissue volumes and regional grey matter changes. 

4. Discussion 

Analysis of subcortical structures revealed significant differences in 
the shape and volume of the brainstem of FMS patients. Voxel based 
morphometry analysis also demonstrated local grey matter reductions 
in the brainstem as well as in the left precuneus of FMS patients. The 
volume reduction of the brainstem of FMS patients showed a significant 
correlation with MTPS scores, indicating that this structure may play an 
important role in the pathogenesis or maintenance of deep tissue pain 
and hyperalgesia in FMS. In addition, total grey matter volume in FMS 
patients also demonstrated a negative correlation with MTPS scores 
which further suggests that grey matter reductions in FMS are relevant 
to symptom severity. 



Table 2 

Local grey matter volume increases in FMS patient group relative to healthy control group. 



Structure 




MNI [mm] 




k 


Z 


T 


Sig. 




X 


y 


z 










Left SI 


-36 


-47 


69 


39 


3.41 


3.83 


0.005 


Right SI 


39 


-39 


59 


35 


3.37 


3.76 


0.005 



Anatomical structure location as defined by the Harvard-Oxford atlas, t-maxima locations 
according with MNI x, y, z coordinates in millimetres. L, left; R, right; k, number of voxels; 
T, peak t values; Z, peak z values; Sig., cluster-level FWE corrected P values; SI, primary 
somatosensory cortex. 
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Fig. 4. Scatter plot of manual tender point scale (MTPS) scores and total grey matter 
volume (voxels) as indicated by VBM8 analyses in the FMS patient group. The linear re- 
gression line is also shown. 



In the FMS patient group a reduction or 'sucking in' of shape was 
evident, located in the left lateral medullary funiculus extending from 
the inferior tip of the medulla to the level of the inferior olivary body. 
Patients also exhibited significant reductions in mean total volume 
of the brainstem. The uncorrected results (Fig. 1) show that similar 
patterns of vertex alterations were seen bilaterally in the FMS patient 
group at the level of the medulla. However, only left sided vertex alter- 
ations were evident following FDR correction in FIRST. Further investiga- 
tion, utilising a larger patient group, is needed to investigate whether 
alterations in the brainstem of FMS patients could actually occur bilater- 
ally. Alternatively, unilateral shape alterations may be sufficient to 
cause widespread symptoms. Small unilateral brainstem infarction has 
been shown to induce major cognitive problems (Garrard et al, 2002; 
Hoffmann and Cases, 2008), and unilateral stimulation of periventricular 
grey can relieve lower back pain (Rasche et al., 2006). Brainstem struc- 
tures likely play an important role in the maintenance of central sensiti- 
sation in FMS (Graven-Nielsen and Arendt-Nielsen, 2010; Petersel et al., 
201 1 ), and it is possible that even unilateral structural alterations could 
be sufficient to contribute to development of widespread central sensiti- 
sation. As no correlation was seen between brainstem vertex positions 
and brainstem volume in the FMS patient group following FDR correc- 
tion, it is not possible to infer that shape alterations in the brainstem of 
FMS patients were directly related to volumetric reductions. However, 
this lack of correlation may be due to the sample size, and this relation- 
ship also warrants further investigation with a larger FMS patient cohort. 

The inward shape alterations evident on the surface of the medulla 
in FMS patients indicate a possible reduction in the nuclei underlying 
this region. Specifically, a reduction in the volume of the reticular forma- 
tion nuclei could lead to the shape alterations seen on the brainstem 
surface. Reticular formation nuclei are involved in many systems 
which demonstrate dysfunction in FMS such as homeostatic regulation, 
postural control and sleep cycle control (Naidich and Duvernoy, 2009). 
This formation also contains the rostral ventromedial medulla which is 
the major relay point between periaqueductal grey and dorsal horn 
neurons, and the pathway of these three structures constitutes the pri- 
mary descending pain modulation pathway in the brain (Lovick, 2008; 
Naidich and Duvernoy, 2009). The nucleus raphe magnus, a group of 
nuclei which also receive afferent input from the periaqueductal grey 
and have efferent connections with dorsal horn neurons, is also located 
in the region, and this structure is also important for descending pain 
modulation. 

Previous studies have demonstrated that descending pain modula- 
tion is impaired in FMS (Julien et al., 2005; Kosek and Hansson, 1997). 



Functional MRI studies of FMS patients undergoing noxious stimulation 
indicate reduced activity in regions of the brainstem associated with de- 
scending pain inhibitory mechanisms (Jensen et al., 2009), and reduced 
connectivity between the brainstem and pain processing regions 
(Jensen et al., 2012). A VBM study also recently reported grey matter 
reductions in the brainstem of chronic pain patients, including FMS 
patients (May, 2009). Functional imaging of the brainstem and spinal 
cord during innocuous and painful touch suggests that dysfunctional 
descending pain modulation may be important in the experience 
of allodynia pain (Ghazni et al., 2010), and such dysfunction may facili- 
tate the development of chronic pain in FMS and other chronic pain 
disorders (Staud, 2011a). In healthy participants, selected regions of 
brainstem showed activation increases during stimulation of a second- 
ary hyperalgesic area of skin (Zambreanu et al., 2005), and a role of 
brainstem areas has also been suggested for anticipatory mechanisms 
involved in the maintenance of chronic pain (Fairhurst et al., 2007). 
Shape alterations and volumetric reductions in the brainstem may 
relate to dysfunctional pain modulation in FMS. However, a recent 
VBM study of chronic fatigue syndrome, a disorder closely related to 
FMS, also suggests structural grey matter reductions in the brainstem 
(Barnden et al., 2011) (although these changes appear to involve dif- 
ferent sub-regions of brainstem). 

In a complex syndrome such as FMS it is difficult to surmise specific 
clinical implications of local abnormalities from a cross-sectional study 
as extraneous factors could also influence structural changes. FMS pa- 
tients do less exercise and exhibit higher levels of obesity than healthy 
people (Mork et al, 2010), and both of these factors were previously 
shown to influence regional and total brain volumes in normal popula- 
tions (Colcombe et al., 2006; Erickson et al, 2009). Duration of medica- 
tion use was also previously shown to correlate with specific structural 
alterations in FMS patients (Burgmer et al, 2009). Therefore, whilst 
every effort was made to recruit a homogenous group of FMS patients, 
it is impossible to ensure that the presence of FMS was the only factor 
influencing brain structure alterations. However, as the shape and 
volume alterations seen in the brainstem of patients did not correlate 
with BDI or FIQ scores, we infer that changes to this structure in FMS pa- 
tients do not relate to changes in mood associated with the disorder. 
Shape changes in the lower brainstem of FMS patients cannot be at- 
tributed to Chiari I malformation which may be associated with FMS, 
or at least with a subgroup of patients (Heffez, 2011; Thimineur 
et al., 2002; Watson et al., 201 1 ). However, we cannot exclude the pos- 
sibility that a combination of multiple anatomical factors in the 
cervico-spinal region might have contributed to the shape change 
seen in lower brainstem of FMS patients in the present study. 

Using VBM analysis, decreased local grey matter volumes were iden- 
tified in the brainstem of FMS patients located in the left ventral aspect 
of the basilar pons. This region does not correspond to the difference 
seen using shape analysis, but contains the pontine nuclei (Naidich 
and Duvernoy, 2009) which are involved in the control and regulation 
of motor activity. As FMS affects dynamic balance control and is associ- 
ated with clumsiness and a high prevalence of falls (Jones et al., 2011), 
grey matter reductions in this region could relate to reported dysfunc- 
tional postural control. The left precuneus also demonstrated local 
grey matter reductions in FMS patients relative to healthy control sub- 
jects. Precuneus is an important structure in endogenous pain modula- 
tion which demonstrates increased connectivity with periaqueductal 
grey and posterior cingulate cortex during experimental pain relief 
(Zyloney et al., 2010), and it is also a known structure of the default 
mode network (Greicius et al., 2003) which can be disrupted by ongoing 
chronic pain (Baliki et al., 2008), and specifically FMS (Napadow et al., 
2010, 2012). Thus, grey matter reductions seen in the precuneus of 
FMS patients may contribute to dysfunctional descending pain modula- 
tion generated in the posteromedial cortex of FMS patients, or to disrup- 
tion of the default mode network. 

FMS patients exhibited bilateral increases in grey matter of somato- 
sensory cortices. Primary somatosensory cortices are important in the 
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discriminative aspects of pain processing (Lee et al, 2008), and aug- 
mented activations were previously seen in somatosensory cortices 
of fibromyalgia patients during noxious pressure stimulation (Gracely 
et al, 2002), and brushing of the forearm (Fallon et al, 2013). 
Neuroplastic alterations in grey matter of chronic pain patients could re- 
late to underlying mechanisms of central sensitisation (Woolf, 2011), 
and dynamic re-organisation of neural networks may occur as tonic 
pain persists in FMS (Staud, 2011b). A longitudinal study of chronic 
pain patients with severe osteoarthritis undergoing hip arthroplasty 
has shown that grey matter alterations normalise post-operatively 
with pain reduction (Gwilym et al., 2010), and volume increases in pri- 
mary somatosensory cortex may ensue in healthy people after just eight 
days of noxious stimulation lasting for only 20 min (Teutsch et al, 
2008). Therefore, grey matter increases in somatosensory regions of pa- 
tients may represent dynamic changes related to impaired sensory pro- 
cessing and discrimination in FMS. However, as experimental pain does 
not accurately reflect the experience of ongoing chronic pain and, due to 
the present uncertainty regarding the precise contribution of peripheral 
afferent input to FMS pain, the cause of these alterations in FMS patients 
requires further study. 

Subcortical morphological alterations in FMS patients are pertinent 
to the disorders pathophysiology and merit further investigation. The 
novel technique of shape analysis used in the present study indicates 
morphological abnormalities in structures relevant to endogenous 
pain modulation and a variety of FMS symptoms. Vertex analysis in 
FIRST is limited to the location of shape changes on the surface of each 
structure, and volume analyses only reveal global volume differences 
for each structure. However, as brainstem volume reductions correlated 
with MTPS scores in the FMS patient group, structural alterations in the 
present study may corroborate previous evidence for dysfunctional de- 
scending inhibition in brainstem of FMS patients (Jensen et al., 2009, 
2012). Therefore, the findings represent a step towards understanding 
the structural underpinning for physiological mechanisms of FMS pain. 

In future, high-resolution imaging studies should prioritise struc- 
tural analysis of cervico-spinal region in FMS patients. Ultra high- 
resolution Tl -weighted MR scanning has already been utilised to iden- 
tify nuclei and white matter tracts in the brainstem of ex vivo humans 
(Soria et al., 2011 ), and future technical improvements may eventually 
allow for in vivo evaluation of brainstem cytoarchitecture in FMS pa- 
tients. Future studies could also further investigate whether subcortical 
shape differences in FMS patients could correlate with functional ab- 
normalities such as decreased functioning of descending nociceptive 
control systems during experimental pain, dysfunctional postural con- 
trol or sleep disturbance. Possible structural and functional alterations 
to the pain modulatory system in FMS may also be affected by current 
treatments (such as pharmacological interventions, self-management 
programmes, or non-invasive brain stimulation techniques), and a 
longitudinal investigation could improve understanding of both the 
pathophysiology of FMS and the mechanisms of such interventions. 
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